Deflection of Rotating Symmetric Molecules by Inhomogeneous Fields 
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We consider deflection of rotating symmetric molecules by inhomogeneous optical and static elec- 
tric fields, compare results with the case of linear molecules, and find new singularities in the distri- 
bution of the scattering angle. Scattering of the prolate/oblate molecules is analyzed in detail, and 
it is shown that the process can be efficiently controlled by means of short and strong femtosecond 
laser pulses. In particular, the angular dispersion of the deflected molecules may be dramatically 
reduced by laser-induced molecular pre-alignment. We first study the problem by using a simple 
classical model, and then find similar results by means of more sophisticated methods, including the 
formalism of adiabatic invariants and direct numerical simulation of the Euler-Lagrange equations 
of motion. The suggested control scheme opens new ways for many applications involving molecular 
focusing, guiding, and trapping by optical and static fields. 

PACS numbers: 33.80.-b, 37.10.Vz, 42.65.Re, 37.20.+j 



I. INTRODUCTION 

Formation of aligned or oriented molecules has long 
been of interest in chemistry and physics. Modern ap- 
plications of aligned and oriented molecules, such as 
high-harmonic generation laser pulse compression 
, nanolithography Q , control of photodissociation and 
photoionization [4[ , and quantum information processing 
0, have motivated the development of all-optical tech- 
niques for aligning molecules under field-free conditions. 
A major advance has been the use of linearly polarized, 
nonresonant ultrashort laser pulses to align molecules by 
an impulsive Raman mechanism Q. Short laser pulses 
excite rotational wave packets, which results in a consid- 
erable transient molecular alignment after the laser pulse 
is over, i.e. at field-free conditions (for recent studies of 
laser induced field-free alignment of non-linear molecules, 
see e.g. These methods also provided tools for 

more general control of the molecular dynamics. In par- 
ticular, we have recently shown that laser induced field- 
free alignment may affect dramatically the scattering of 
linear molecules from inhomogeneous optical and static 
electrical fields This opens new ways for many 

applications in molecular optics, involving molecular fo- 
cusing, guiding and trapping by optical and static fields 

As most molecules are not linear, a more general the- 
ory is required for describing scattering of symmetric and 
asymmetric molecules by inhomogeneous fields. In this 
work, we focus on the scattering of symmetric molecules 
by optical and static electrical fields, and investigate dis- 
tinctive features in the distribution of scattering angles 
due to molecular rotation. In particular, we investigate 
rotational rainbows in scattering, and find a new one in 
addition to those described in our previous study on the 
scattering of linear molecules [ll|. The difference be- 
tween the symmetric and the linear molecules can be eas- 
ily seen at the classical level. Having three non- vanishing 
moments of inertia, symmetric molecules rotate in a more 
complicated fashion [17], and unlike the linear molecules, 
their rotation is not constrained to a plane perpendicular 



to the angular momentum. We will show that the scatter- 
ing of symmetric molecules can be significantly affected 
and controlled by preshaping molecular angular distribu- 
tion before the molecules enter the interaction zone. This 
can be done with the help of numerous recent techniques 
for laser molecular alignment, which use single or mul- 
tiple short laser pulses (transform limited or shaped) to 
align molecular axes along certain directions. This paper 
is based on classical considerations, since we have demon- 
strated in the past a good correspondence between the 
quantum mechanical and classical treatments [Til Il2j of 
the strong-field rotational control. The first part of the 
paper, i.e. Sees. Hill VII are devoted to the molecular de- 
flection by optical fields. The last part of the paper, i.e. 
Sec. IVIII describes the deflection of symmetric molecules 
by a static electric field. In Sec. [TT| we introduce the 
general description of the molecular deflection by opti- 
cal fields. In Sec. IIII1 a simple model of a symmet- 
ric molecule is presented, for which rotational rainbows 
are predicted in the distribution of the scattering angle. 
In Sec. IIVI we generalize the simplified model to ther- 
mal conditions. In Sec. |Vj we consider examples of real 
molecules, and compare between the oblate and prolate 
cases. In Sec. IVI1 laser-induced field-free alignment (for 
various polarizations) is implemented in order to control 
the scattering process. We expand the study to deflec- 
tion by a static electric field in Sec. IVIII and finally, we 
summarize in Sec. IVIIII 



II. DEFLECTION OF MOLECULES 

Although our arguments are rather general, we follow 
for certainty a deflection scheme similar to one of the 
experiment by Stapelfeldt et al [l4[ who used a strong 
IR laser to deflect beams of linear molecules, and then 
addressed a portion of the deflected molecules (at a prese- 
lected place and time) by an additional short and narrow 
ionizing pulse. Consider deflection (in Z direction) of 
a molecule moving in X direction with velocity vx and 
interacting with a focused nonresonant laser beam that 
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propagates along the Y axis. 

The spatial profile of the laser electric field in the XZ- 
plane is: 

E = E cxp[-(X 2 + Z 2 )/a;2]exp[-21n2t 2 /^ 2 ]- (1) 

We consider a laser field polarized in the Z direction, 
and find the interaction potential of the laser with the 
asymmetric molecules to be 18]: 

U a = -jE 2 (a ab cos 2 9 + a cb sin 2 6 sin 2 x + a b ) , (2) 

where E is defined in Eq. [TJ and 9 and x are the Eulcr 
angles. a a ,b,c are the polarizability components of the 
molecule, along its principal axes (a, b, c), and a 13 = at — 
oij. Here, the Euler angles, as well as the principal axes of 
the molecule are defined as in Fig. 2 of Ref. [19j. In the 
case of symmetric molecules, we can define the parallel 
and perpendicular polarizability components (ay and a±, 
respectively) as: a a = an and cvf, — a c = au_- Using Eq. 
[2l one obtains the interaction term: 

U s = -\e 2 [(a|| - a ± ) cos 2 9 + a±] , (3) 

which has the same form as for the linear molecules. 

A molecule initially moving along the X direction ac- 
quires a velocity component vz along Z-direction. We 
consider the perturbation regime corresponding to a 
small deflection angle, 7 « Vz/vx- We treat Z as a 
fixed impact parameter, and substitute X = Vxt. By 
doing this, we concentrate on the molecules reaching the 
focal spot at the moment of the maximum of the deflect- 
ing pulse, like in Refs. [13, [Hj]. The deflection velocity 
is given by: 

VZ = ^Il Fzdt = -^Il^ Us ) Z dt - (4) 

Here M is the mass of the molecules, and Fz is the de- 
flecting force. The time-dependence of the force Fz (and 
potential U s ) in Eq.Q comes from three sources: pulse 
envelope, projectile motion of the molecule through the 
laser focal area, and time variation of the angle 9 due to 
molecular rotation. For simplicity, we start with the case 
of a relatively weak deflecting field that does not affect 
significantly the rotational motion. Such approximation 
is justified, say for Benzene molecules with the rotational 
temperature T = 5K, which are subject to the deflecting 
field of 3 • 1{) 9 W / cm 2 . The corresponding alignment po- 
tential U w — j (aii — E 2 w 0.02 meV is an order of 
magnitude smaller than the thermal energy fc^T, where 
ks is Boltzmann's constant. This assumption is even 
more valid if the molecules were additionally subject to 
the aligning pulses prior to deflection. 

Since the rotational time scale is the shortest one in 
the problem, we average the force over the fast rotation, 



and arrive at the following expression for the deflection 
angle, 7 = v z /v X - 

7 = 7o [(ay - ax)A+ ol] • (5) 



Here A = cos 2 8 denotes the time-averaged value of 
cos 2 9. This quantity is different for different molecules 
of the incident ensemble, which leads to the randomiza- 
tion of the deflection process. The constant 70 presents 
the average deflection angle for an isotropic molecular 
ensemble: 

f-4Z\ 
70 AMv 2 x {coo J 

( 2w 2 ln2V 1/2 / 2Z 2 \ , . 



III. SIMPLE MODEL 

We consider a symmetric rotor, as in Fig. [1] (the co- 
ordinates are defined in the figure). This rotor precesses 
about its own axis, and the axis itself rotates about the 
angular momentum J (the angle 6i between the molecu- 
lar axis and J remains constant). As seen from Eq. [31 the 
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FIG. 1: The symmetric molecule model. X,Y,Z are the lab- 
oratory coordinates, and 9j is the angle between the angular 
momentum J and the Z axis. The x,y,z coordinates are 
obtained by rotating the X, Y, Z coordinates about Y by an 
angle 6,j, such that the z\\J . The molecule precesses about 
its own axis, and the axis itself rotates in a conical trajectory 
about J. 

time-averaged interaction of the symmetric molecule and 
the laser field is linearly related with A (A = cos 2 9, as 
was defined before). In what follows, we will estimate the 
magnitude of A. In the x, y, z coordinates, the molecular 
axis direction is given by: 

r x ,y,z = {sin 9 i cos Lot, sin 9 ism tot, cos 9 i), (7) 
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where uj is the frequency of rotation. The relation be- 
tween the x, y, z and the X, Y, Z coordinates is: 



cos 8 j —sin 
1 
sin 9 j cos 9j 




(8) 



Therefore, the molecular axis in the X, Y, Z coordinates 



is: 



cos 9 j — sin Oj 
rx.Y.Z = j 1 | 
sin Oj cos 9 j 



x,v,z 



(9) 



and the projection of the molecular axis on the Z axis is 
given by: 

?xy z ■ Z = cos 9 = sin 9 j sin 9i cos ut + cos 9j cos 9i , 

(10) 

where 9 is the angle between the molecular axis and the 
Z axis. Taking the square, and averaging it over time, 
we obtain: 



A — „™2 a _„;„2 a i „™2fl_„„„2, 



■ sm t)j sm 



cos uj cos 



(11) 



To find singularities in the distribution of scattering an- 
gles, we determine the extremal points ^ = 0, and 
find that: 



OA 

09, 



i.j 



sin 29i = and sin 26 j = 

or (12) 



sin 0j = | and sin #j 



According to the first line of Eq. [121 the extremal an- 
gles 9i and 6*j can be either 0, tt/2, it. Considering smooth 
distributions of the angles 9i,0j, the regions near 0, tt give 
a negligible contribution to the distribution of A because 
of the vanishing phase space volume. Therefore we con- 
clude that the distribution of A may have singularities 
at: 

A(6i = it/2, 9j = it/2) = 1/2 

.4(sin 2 6i = 2/3, sin 2 9j = 2/3) = 1/3. (13) 



Using this equation, and also the expression for the total 
energy defined by Eq. Q3J one can easily find the angles 
9i , 9j introduced before: 



sinf 



COS Oj 



J 



2I C E 



(16) 



For a thermal molecular ensemble with temperature 



T, we define dimensionless momenta Pi 



- P v cos 6 



pi = 
X ^I a K B T 



and P' 



and calculate the dis- 



tribution of the time-averaged alignment factor A as: 



d9(0)sm9{0)d(j>(0)dx(0) 



x dPl){0)dP^{0)dP^0)6{A - cos 2 9) 

x F(0(O),0(O),P^(O),P^(O),P;(O)), (17) 



where 



F 



Qr 



■ exp 



-L(pZ + p? + p£) 



(18) 



and Q ro t is the partition function. Here cos 2 9 is given 
by Eq. [TTJ 

A Monte-Carlo simulation was performed, in order to 
account for the random initial orientation and momenta 
of the molecules in the thermal ensemble. The results 
are shown in Fig. [21 where we see singularities at the 
points A = 1/3, 1/2, in full agreement with the previous 
speculations. 




IV. SCATTERING OF A THERMAL 
ENSEMBLE 

Consider a symmetric rigid rotor described by a La- 



grangian: 



2I r 



2I r sin 



2 , 



p2 

w.: < 14) 



where Ii, = I c , and I a is the moment of inertia about 
the molecular symmetry axis. Pg, P^ and P x are the 
canonical momenta. The angular momentum J can be 
found by: 



J 2 =P 2 



P 



{P 4> -P x cos9f 
sin 2 9 



(15) 



FIG. 2: Distribution of A for a symmetric molecule. T — 5K, 
and the following values of the moments of inertia were used: 



h = i-io- 



2-10" 4B fcc/-m / and I c = 2-10" 



m . The calculations are based on Eq. 1111 and a Monte-Carlo 
simulation. 

We have also verified the results presented in Fig. [Hby 
using two other methods. First, based on the separation 
of the rotational and translational time scales, we used 
adiabatic invariant technique [13, H3, HI! to calculate A 
for a symmetric molecule moving in inhomogeneous field, 
and utilized these data for calculating the deflection an- 
gles. In the second method, we numerically solved the 
Euler-Lagrange equation of motion, calculated cos 2 9(t), 
and averaged it over a long enough time to reach con- 
vergence [12[- In both methods, the final results were 
similar to those shown in Fig. [21 
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V. REAL MOLECULES: OBLATE VS. PROLATE 

We are ready now to apply the developed procedure 
to real molecules. We start with Benzene, i.e. an oblate 
molecule. For this molecule, the symmetry axis is a, and 
I a = 2-I b , where I b = I c = 1.474- 10~ 45 fcg-m 2 (moments 
of inertia and polarizability data are taken from [22|). 
The distribution of A for a Benzene molecule is given in 
Fig. [31 It is seen that the distribution peak near .4=1/3 
is the dominant one in the case of Benzene. 




0.2 0.4 0.6 0.8 1 
A 



FIG. 3: Distribution of A, for the Benzene molecule. The 
calculations are based on Eq. 1111 

As another example, we consider Methyl Iodide 
(CH3I), which is a prolate symmetric molecule. Again, 
the symmetry axis is a, but the moments of inertia are: 
I a = 5.4071 ■ 10- 47 kg • m 2 , I b =I c = 1.118 • lO'^kg ■ m 2 
(moments of inertia data are taken from [23j , and polar- 
izability data are taken from [24j). The corresponding 
distribution is presented in Fig. [4] For this molecule, 
the most visible peak in the distribution is around the 
A = 1/2 point (quite similar to the case of linear 
molecules |11|). 
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°0 0.2 0.4 0.6 0.8 1 
A 

FIG. 4: Distribution of A, for the Methyl Iodide (CH 3 I) 
molecule. The calculations are based on Eq. 1111 



VI. CONTROLLING THE SCATTERING OF 
SYMMETRIC MOLECULES BY A LASER FIELD 

Here we consider the laser deflection of symmetric 
molecules, that are prealigned by a femtosecond laser 
pulse before entering the deflecting field. The polariza- 
tion of the prealignment laser pulse is not necessarily in 
the direction of the deflecting field, as is shown in Fig. 
[5) After the end of the prealigning pulse, the quantities 



Z 




FIG. 5: The deflecting field is in the laboratory Z direction, 
and the prealignment pulse (denoted as e) is tilted by an angle 
a with respect to it. 



Pe(0) and P^(0) are replaced by: 

P 8 (0) -> P e (0) + P S (sin 2 a sin 2(9 cos 2 

+ sin 2a cos 29 cos <fi — sin 29 cos 2 a) 
P (O) -> P (O) - P s (sin 2 a sin 2 sin 20 

+ isin2asin26>sin(/>), (19) 

where P s = Ph and P = (l/4h)(a\\ - a±) J e 2 {t)dt. 
The dimensionless parameter P is the so called "kick 
strength" that equals to the typical amount of angular 
momentum, in the units of ft, supplied by the pulse to 
the molecule. For Benzene, Aa = ay — a± is a negative 
quantity, therefore P is negative. 

We analyzed the effect of prealigning pulses on the 
scattering of two molecules considered in the previous 
section. The kick strength was taken as large as |P| = 25. 
This corresponds to a prealignment pulse with a peak 
intensity of 4.13 • l0 12 W/cm 2 and a duration of 0.5ps 
(FWHM) for Benzene. For CH 3 I, P = 25 corresponds 
to the pulse intensity of 2.41 • 10 14 W/cm 2 and a duration 
of 20 fs (FWHM) . Fig. ©shows distribution of A both for 
Benzene and CH3I , after the molecules were prealigned 
by a laser pulse tilted at different angles a (0° , 45° and 
90°) with respect to the deflecting field. In the prolate 
molecule case, the distribution is somehow similar to the 
distribution of a linear rotor that was kicked by the pre- 
alignment pulse. When the prealignment polarization is 
parallel to the deflecting field (a = 0°), the vector of the 
angular momentum of the molecules is preferentially con- 
fined to the XY plane after the pulse, i.e. 9i sa 9j tt/2, 
which corresponds to A — 1/2. In this way, the molecules 
experience the maximally possible time-averaged deflect- 
ing force which is the same for all the particles of the 
ensemble. As the result, the dispersion of the scattering 
angles is reduced dramatically (Fig. CH3I, a = 0) [PH ]. 
For other polarization angles, two peaks are observed i.e. 
at A = 1/2 and (cos 2 a)/2 (Fig. H CH 3 I,a = 45,90) 
[l3l |. For oblate molecules, the difference from the linear 
molecules is more emphasized due to the dominant peak 
at A = 1/3. 
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FIG. 6: Benzene ((a)-(c))and CH 3 I ((d)-(f)) A distributions, 
after the ensembles were prealigned by the means of a fs laser 
pulse, at different angles a ((a) and (d) correspond to a — 0°; 
(b) and (e) correspond to a — 45°; (c) and (f) correspond 
to a = 90°). T = 5K, P — -25,25 for Benzene and CH 3 I, 
respectively. The right column (Prolate) is very similar to 
the linear molecule, where the left column (Oblate) shows the 
peak at 1/3, in addition. 



VII. DEFLECTION OF SYMMETRIC 
MOLECULES BY STATIC ELECTRIC FIELDS 

In this section we expand the research to scattering 
of symmetric molecules by static inhomogeneous electric 
fields. We follow here the same line of thought presented 
in our recently published paper [l2| , where we considered 
a deflection by a static electric field of a beam of linear 
molecules. Briefly, in the case of deflection by a static 
electric field, the field interacts not only with the molec- 
ular polarizability, but also with its permanent dipole 
moment. Thus, Eq. [3] is modified as: 

U s = ~^E 2 [(ay - aa) cos 2 6> + aa] -pEcosO. (20) 

According to Sec. [TTl the important variables defining 
the deflection are given by the time-averaged values of 



cos9 and cos 2 9, i.e. 



cos 9 



A 1 

A 2 = cos 2 9. 



(21) 



As in [12j , we estimate the strength of the deflecting field 
by means of the following parameters: C = E 2 Aa/ (kgT) 
and D = fj,E/(k B T). For instance, for a 10 K CF 3 H 
molecular beam that is deflected by a static field of 
2 x lOV/m, the corresponding parameter values are 
C = -9.7 x 1CT 5 and D = 0.8 (the typical electrical 
properties of the CF^H molecule were taken from 25]). 
In this case, the interaction with the induced polariza- 
tion may be neglected. We consider only the interaction 



with the permanent dipole moment of the molecule, and 
calculate A\ and A2 by means of the adiabatic invari- 
ants method (the derivation is given in [20]). The re- 
sults are presented in Fig. [7J The distribution of the 
deflection angle is linearly related with the distribution 
of A\ (Fig. [7^,), however we provide also the distribution 
of A2, both for educational purposes and for practical 
ones, since one can think of an experiment to measure 
the Ai distribution, i.e. by means of an additional opti- 
cal field [l2j]. Fig. [7J is very similar to Fig. 5 from [l^l (in 




FIG. 7: Statistical distributions of (a) .Ai and (b) A2 for a 
thermal beam of CF 3 H molecules in a moderate deflecting 
field (2 x lOV/m). T = 10K, C = -9.7 x 1(T 5 and D = 0.8. 

both cases the D's values are of the same order of magni- 
tude). Therefore, the interpretation of the results given 
in [l2j is applicable to the present case too. Briefly, the 
tail of the distribution is formed by relatively low-energy 
molecules that are angularly trapped by the deflecting 
field (Ai and Ai are close to 1). Some of the molecules 
in the ensemble are rotating in vertical planes that con- 
tain the deflecting field polarization. They contribute to 
the negative peak of Fig. [7^,, and the rainbow peak in 
Fig. [7b around the value Ai = 0.28 (for slightly trapped 
molecules) and A2 = 0.5 (for free molecules). The main 
difference between Fig. [7p and Fig. 5b from 12] is the 
suppressed peak value at A2 = 0.5 in Fig. [7b, as well as 
the less clear peak at A2 — 0.28. The suppression of these 
peaks is due to the important difference between linear 
and symmetric molecules, namely the possibility of the 
symmetric molecules to rotate around their own symme- 
try axis. Therefore, unlike the case of linear molecules, 
which rotate both with and against the electric field, a 
portion of the symmetric molecules in the ensemble may 
interact only along the symmetric field (and hence be 
more deflected by it). We also examine the case of a 
stronger deflection, i.e. by means of a 2 x 10 8 V/m de- 
flecting field, where C = -9.7 x 10" 3 and D = 8. The 
results are presented in Fig. [51 where a larger portion 
of the molecules is trapped by the electric field. This 
Fig , may be compared to the quite similar Fig. 6 from 
[l2j. Finally, we apply a prealigning laser pulse (P = 25) 
polarized parallel to the direction of the deflecting field 
(i.e. in the Z direction). The results are given in Fig. [9] 
and are very similar to the results presented in Fig. 8 of 
[l2^ . Due to the initial prealignment, the molecules gain 
a strong rotational energy and the interaction with the 
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FIG. 8: Statistical distributions of (a) A\ and (b) Ai for a 
thermal beam of CF3H molecules in a strong deflecting field 
(2 x 10 8 V/m). T = 10K, C = -9.7 x lCT 3 and D = 8. 

deflecting field is almost cancelled out (Fig. |^). Further- 
more, the rainbow at A2 = 0.5 is enhanced due to the 
strongly rotating molecules in the vertical planes. Thus, 




0.2 0.4 0.6 0.8 1 



FIG. 9: Statistical distributions of (a) Ai and (b) A2 for a 
prealigned (P = 25) thermal beam of CF3H molecules in a 
strong deflecting field (2 x W 8 V/m). T = 10A", C = -9.7 x 
KT 3 and D = 8. 

we show that by means of laser- induced prealignment, a 
dramatic control of the scattering distribution is possi- 
ble. In particular, the interaction between the molecule 
and the static field can be switched-off (and on). 



new singularities in the distribution of the scattering 
angle. Depending on the shape of the molecule, i.e. 
oblate/prolate geometry, the singularities were shown to 
be suppressed or enhanced compared to the case of linear 
molecules. We developed a simple model for the scatter- 
ing process, and found a good correspondence between 
its predictions and the ones of more complicated treat- 
ments, including adiabatic invariants formalism and di- 
rect simulation of the Euler-Lagrange equations of mo- 
tion. Our results indicate that laser-induced prealign- 
ment provides an effective tool for controlling the deflec- 
tion of the molecules, and it may be used for increasing 
the brightness of the scattered molecular beam. Fur- 
thermore the interaction between the molecules and a 
static electric field can be switched-off (or switched-on) 
by means of the prealignment. This may be important for 
nanofabrication schemes based on the molecular optics 
approach [l8j]. Moreover, molecular deflection by nonres- 
onant optical dipole force is considered a promising route 
to separation of molecular mixtures (for a recent review, 
see Ref. (HI). Narrowing the distribution of the scat- 
tering angles may substantially increase the efficiency of 
separation of multicomponents beams, especially when 
the prealignment is applied selectively to certain molec- 
ular species, such as isotopes [13] or nuclear spin iso- 
mers 0, [HI ■ The predictions of our theory may en- 
courage others to design more complicated schemes for 
controlling the process, e.g. by using several prealign- 
ment pulses with variable polarizations. Future studies 
may also focus on the scattering of asymmetric molecules 
having more complicated deflection potential. Moreover, 
the same mechanisms may prove efficient for controlling 
inelastic molecular scattering off metalic and dielectric 
surfaces. These and other aspects of the present problem 
are subjects of an ongoing investigation. 
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